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ULTRARAPIDLY DECREASING ULTRADIFFERENTIABLE FUNCTIONS, 
WIGNER DISTRIBUTIONS AND DENSITY MATRICES 

JEAN-MARIE AUBRY 
Abstract 

Spaces <Soj , <S{jj } i ) of ultradecreasing ultradiffcrentiable (or for short, ultra-5) functions, depending on a 
weight e"' 1 ', are introduced in the context of quantum statistics. The corresponding coefficient spaces in the 
Fock basis are identified, and it is shown that the Hermite expansion is a tame isomorphism between these 
spaces. These results are used to link decrease properties of density matrices to corresponding properties of the 
Wigner distribution. 

1. Introduction 

1.1. Quantum states 

A quantum state is a vector of unit norm in some Hilbert space, for instance L 2 (M). In 
quantum optics (our subjacent model throughout this paper), the simplest states decrease like 
the Gaussian, and have the same regularity (their Fourier transform decreases like the Gaussian 
as well). One can however easily encounter super Gaussian states whose wave function decreases 
still rapidly, but not as fast as the Gaussian, like for instance \f(x)\ < Ce 2 for < (3 < 2, 
and whose Fourier transform decreases at the same rate (the restriction (3 < 2 comes from the 
uncertainty principle, detailed in § I2.2|) . Rapidly decreasing self- Fourier functions, of importance 
in optics [3 [9l 110) , fall under this category. 

The class of compactly supported functions / that have a regularity comparable to — 

lil 13 

Ce ~ is called a class of ultradiffcrentiable functions (in the sense of Bcurling, see [5]). The 
classes of functions that we introduce in §[2] are related but strictly larger, because we replace 
compact support with an "ultrafast" decrease condition. 

When both conditions (|/(x)| < Ce 2 and /(£) < Ce 2 ) arc satisfied, / can be called 
an "ultrafast decreasing ultradiffcrentiable" function. More generally, following Grochenig and 
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Zimmerman |12] , we can replace x i— > ^jj- by a weight function u> satisfying certain conditions 
(Definition [l}. The spaces thus obtained can be viewed as the ultradifferentiable version of the 
Schwartz class S; for short we call them spaces of ultra-S functions. 

1.2. The Fock basis 

The natural basis of L 2 (K) for the problems of quantum optics is the one formed by the 
Hcrmite functions 



where H m (x) '■= ( — l) m e x ~ -§^k£~ x is the m-th Hermite polynomial. The normalization above 
ensures that (h m ,h n ) = 5 m ,n- 

In terms of time/frequency, the Hermite functions arc perfectly balanced, because they are 
eigenvectors of the Fourier transfornQ 



This basis is also called the Fock basis and the h m are the Fock states (the m-photons state 
in the case of light). Note that (Q]) characterizes "self- Fourier" functions as the ones having 
non-zero Hermite coefficients only for m £ 4N. 

A few more properties will be needed. For instance, h n is also an eigenvector of the Hermite 
operator 



Some decrease properties of h n (x) are stated in Lemma 15771 

1.3. Outline of the results 

Despite their simple definition, the study of the functional properties of Hermite series 
is in many cases a difficult problem. The L p case, for instance, is the object of a book by 
Thangavelu [26] . A question similar to ours has also been investigated numerically by Boyd 
in [4]. The case u>(x) = Cx was treated by Janssen and van Eijndhoven [17j . More anecdotically, 
the decrease of the Hermite coefficients of a certain function has been related to properties of 
the Riemann £ function by Grawe in Finally, Langenbruch |21j has treated the problem 
starting from a slightly different definition for the function spaces (conditions of the type 
||a; Q £)^/|| < M a+ p); in fact a part of his proofs could be adapted and reused here. 




{-i) m h. 



(1) 



- lh) h n(x) = (2n + l)h n {x). 



tThe definition that wc use is /(£) := {2tx)~^ J f(x)e-^ x dx 
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Our first goal in this paper is to characterize the class of ultra-5 functions introduced above 

in terms of Hcrmite coefficients; this is done, up to a multiplicative constant on ui, in §[3J The 

main result in this paper is Theorem [T] in § 13.11 

The next question occurs when one considers a so-called mixed state: that is, the statistical 

result of a physical experiment where there is an uncertainty on the quantum state itself. Let us 

say that, with probabilities po, pi, . . . , the experiment produces orthogonal states ipo, ipi, ■ ■ ■ ■ 

All this information is contained in the operator 

! 

where (.,ipi) i^i is the projector on the subspace generated by ipi. If more than one of the Pi 
is non-zero, then the state is called mixed. In other words, a pure state is one that can be 
described by a single wave functions: for instance ijj = y / p7?/'i- 

It is easily seen, as in Leonhardt's reference book [22] , that p is Hermitian, that its eigenvalues 
are nonnegative with a sum tr(p) = 1, and that these conditions are sufficient to make a mixed 
quantum state. 

For a more graphic representation, one can also consider the Wigner distribution associated 
to p. It is a function of two variables 5,p£R (or rather, a function of q + ip G C), defined as 

\- P 



® P (q,p) - 



lpa W<Z-f)<M<7+f)^. (2) 



For some reasons, it is sometimes easier to work with 



(JL JL\ 



2 P \V2'V2J 

The transform p i— > $ p is an isometry between £(L 2 (R)) and L 2 (C); its inverse is called the 
Weyl transfornx}. The Wigner distribution is very useful for the interpretation of experiments, 
such as in quantum homodyne tomography (op. cit., see also the lectures notes edited by Paris 
and Rehacek |24j). 

For a pure state ip, or in terms of operators, p^ = (.,"0) the Wigner representation of p^ 
is given by a simpler form of ([2]) called the Wigner transform of ip, which is the quadratic form 

1 

2^ 



ip(q-^Mq+^)dx. (3) 



Then there is no difficulty to see that J §(ij},ip)(q,p)dp = \ip(q)\ 2 and that § &(i[i,ip)(q,p)dq = 
h/>(p) ! this is what makes this representation useful in quantum mechanics, for it allows to 



^See 1261 for more insight into this deep connexion with the Heisenberg group representations. 
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recover the probability distribution of an observable and its dual (e.g. positition and momen- 
tum, electric and magnetic field. . . ) as marginal distributions of the Wigner transform of "0 
(even though $(?/>, ip), being in general non-positive, is not a probability density stricto sensu). 

Another consequence is that, generally speaking, a decrease condition on <fr(ip,ip) implies a 
decrease condition on both tp and ip. From the results of §0 this implies a decrease condition 
on the Hermite coefficients of the object of § [4] is to show that the converse is true. 

This result can be extended to mixed states in the following manner: let [p m ,n] be the 
(infinite) density matrix of p in the Fock basis. For instance, in the case of a pure state p = p.,p, 



Pm,n = ("0, h m ) (ip, h n ). (4) 

Assuming that p m , n decreases ultrarapidly when m + n — ► oo, we shall show that <& p also 
decreases ultrarapidly when \p\ + \q\ — > oo. However, in § 14. 21 we see that there are (necessarily 
mixed) states for which $ p decreases ultrarapidly but not p m ,n- 

Another object of interest is 4> p , the 2-dimensional Fourier transform of $ p . In radar tech- 
nology, this function is known as the ambiguity function. Its decrease rate rules the regularity 
of <!>p, which is of importance in statistical estimation, see Butucea et al. [6]. We study the 
decrease and regularity properties of <1> P in § 14.31 

This leads to studying the matrix elements of the Weyl transform: combining j2|) and ([3]) we 
see that, writing 



®m.n(q,p) ■= ®(h m , h n )(q,p) = — 

Ztt 



e ipx h m (q-^)h n (q+^)dx, (5) 



we have just $ p = J2 m .n Pm,n$m,n- Similarly, we define 

$ m Ai,p) : =5*-« (^'^)- 

These functions, named special Hermite functions by Strichartz [26] , are remarkable. Although 
they are not exactly 2-dimensional Hermite functions (the latter being defined by tensor 
product h u ^ v (q,p) := h u (q)h v (p)), special Hermite functions are a special linear combination 
(for u + v = m + n) of those. Consequently, they are also eigenfunctions for the 2-dimensional 
Fourier transform: 

|Q = (_i)^+»$ ro)n . ( 6 ) 



Let us now expose our frame of work. 



ULTRARAPIDLY DECREASING ULTRADIFFERENTIABLE FUNCTIONS 



■3 



2. Functional setting 



2.1. Ultra-S spaces 



Because of ([T]), any function space for which the Fock states (Hermite functions) form an 
unconditional basis must be invariant under Fourier transform. Such spaces can be defined by 
weighted L°° conditions on both / and /, as follows. 



Definition 1. A continuous increasing function uj : [0,oo) 
function if it satisfies 

(i) f2 : t h- > u;(e { ) is convex 

(ii) log(i) = o(w(t)) 

(iii) uj{2t) = 0(w(t)) 



[0, oo) is called a weight 



(iv) limsup 



;(*) 



This definition should be compared with that of Braun, Meise and Taylor [5]. The only 
difference is that their (/?) ^j^-dt < oo is replaced by (iv) which is weaker; this is because 



we relax the compact support condition. As we shall see below, (iv) is related to the uncertainty 



principle; for technical reasons, we shall sometimes need instead the strict inequality 



(v) limsup 



u,(t) 



< 



Given a weight function u> satisfying (i) (iv) we define 



:= {/ G L 2 (R)J(x) = 0(e-*^l>) and /(£) = 0(e~^)} 
Clearly is a Banach space, if equipped with the norm 



ll/IL : = H*)e 



KM) 



4\m 



Remark that A i— » Sx^, is strictly decreasing (in the sense of inclusion) on (0, +oo). Never- 
theless, it will be convenient to state our results in terms of spaces that depend only on the 
rate of growth of u>, that is, on the class ZJ := {Xoj, A > 0}. For this purpose, we introduce the 
ultra-5 spaces 

*%7} : = [J SeuJ 



e>0 



and 



AT<oo 
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Note that the constant in (iii) depends on u only, we shall often use this fact. Also remark that 
in the above union and intersection, the parameter can be taken in a countable set without 
changing anything to the definition. 



2.2. Uncertainty principle 

Bearing the name of Heisenberg, the principle that says that one cannot know with preci- 
sion both the position and the momentum of a particle translates quantitatively in terms of 
simultaneous localization of a function and its Fourier transform; for a comprehensive study 
follow Havin and Joricke |15) . In our setting, this principle says that if w increases too fast, 
then S u is trivial (reduced to {0}). 

More precisely, let us quote the one-dimensional version of a theorem of Bonami, Demange 
and Jaming [3] , in the descent of Hardy |14j , Beurling and Hormandcr 



Theorem (BDJ). Let f € L 2 (R) and N > 0. Then 

\f(x)\\f(0\eW 



(i + M + KI) 



JV 



dxd£ < oo 



if and only if f(x) — p(x)e ^ } where p is a polynom of degree < ^p-- 



This theorem means that if 



:= inf I N, 



— — dxd£ < oo 



R+; 



then the dimension of S u is 



N,. : 



When ^fti e N, then one has to check whether JJ e jjj^j^nz — dxd£, is finite or not, to 



see whether dim(5 w ) 



1 or N " 



2 2 

condition for Su, to be non trivial is that 



(i+*+?) J 

JK 1 XJK. ' 

- respectively. In particular, the necessary and sufficient 



dxdt; 



oo. 



R+ xR+ 

The "critical case" oj{x) = \ is however a fuzzy frontier, because it is actually the limit 



inferior of 



,(x) 



that counts when determining whether the above integral diverges or not. 



Indeed it is easy to construct a weight function u> such that lim inf < | < lim sup (the 



j(x) 



only thing to take care of is (i) a piecewise affine function f2 with increasing slopes will do). 
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This shows that (iv) a fortiori (v) is indeed a restriction to the range of our rcsultsQ There 
might be some room for improvement here. 

Another consequence of Theorem (BDJ) is that a sufficient condition for S^y to be non- 
trivial is given by limsup ^jQ- < oo; for it is that limsup ^W- = 0. 

For comparison, we recall that the fastest possible decrease for the Fourier transform of 
compactly supported functions is of order e _c ' 4 ', see [UH] for instance. 

2.3. Sequence spaces 

The spaces of Hermite coefficients for functions in and Sfa) will be identified in § [3] to 
the following spaces. Let w be a weight function as in Definition [TJ Then 

A w := {(«„) G C N ,a„ = 0(e"^))} 

endowed with the norm ||(a„)|| w := sup„ |a re |e £tf ( v '") is a Banach space, and as previously we 
define 



as well as 



A{„} :— [J A ew 
A<a) ■■= fl A 



Nu 
N<oo 

2.4. Topology 

As a decreasing intersection of Banach spaces, endowed with the projective topology 
(the coarsest topology that makes every embedding Sfc>\ — > Snu continuous) is a Frechet space. 
The same holds for A^j. 

The picture is a little more complicated for Sip\ and A^j . As an increasing union of Banach 
spaces, (A Qa) C Ap u if /3 < a), A{zj} is naturally endowed with the inductive limit topology 
(the finest locally convex topology such that every embedding A euJ — > A/^j is continuous). 
In the standard terminology, A/^n is called a (LB)-space. The same goes for S^y It is a 
classical problem in functional analysis to study the properties of such spaces, in particular 
their completeness. Note that this inductive limit is not strict (in the sense of Kothc [201 
§ 19.4]) because the topology on K aul is not induced by that of Ap u . 

Proposition 2.1. A/^j is a Silva space, thus complete. 



tin fact |(iv)] is only used in Lemma 13.41 to guarantee that the constant K is universal, but any fixed number 
could replace i. Theorem [T] still holds if we replace [(iv)| by a finite limit superior. However, when [(v)] is asked 
for in Propositions l3.2l 14.11 l4~2l then ~ is indeed the critical value. 
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Proof. Recall, as in [8], that a Silva space is a locally convex inductive countable union 
of increasing Banach spaces, such that each embedding map is compact (it is in particular 
a (DFS)-spacc). This is the case for A{^}, as we can see that the union may be taken over 
e { m> 171 e ^0' ^ na ^ ^— ^ " ! an d that the unit ball of A<a. is relatively compact in A * i 



because e ™ = o( e "+ 1 J as n 



Corollary 2.2. S^y is complete. 



□ 



Proof. Anticipating a little, by Theorem [TJ is the image of A/m by a topological 

isomorphism. □ 

2.5. Tame isomorphisms 

The notion of tameness was introduced in [13] . in connexion with the Nash-Moser inverse 
function theorem. It is a natural notion of regularity for linear operators between Frechet spaces 
or inductive limits of Banach spaces. 

Definition 2. Let (Ej,\ \j)jeN and (Fj,\\ be two increasing families of Banach 

spaces, and let E := [jj Ej, F := (J. Fj be the corresponding (LB)-spaces endowed with the 
inductive topologies. A linear mapping T : E — > F is called tame if there are jo G N, C < oo 
such that for all j > jo, there exists C < oo satisfying 

im/)ita < ci/i,.. (7) 

Let (Ej,\ L)jeN and || ||j)jeN be two decreasing families of Banach spaces (or more 
generally, spaces with increasing semi-norms), and let E := |"| . Ej, F : = |"| . be the corre- 
sponding Frechet spaces endowed with the projective topologies. In that case, a linear mapping 
T : E — > F is said to be tame if there are jo G N, C < oo such that for all j > jo, there exists 
C < oo satisfying 

im/)iij < ci/b- (8) 



A linear mapping T is called a tame isomorphism if it is bijective and if both T and T 1 are 
tame. 
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3. Hermite expansion 



3.1. Main results 

In our setting, Q of Definition [2] will apply to spaces of type {uj} and ([5]) will apply to 
spaces of type (u;). 



Theorem 1. Let H : f i-> ((/, /i„))„ 6N . 

If lim sup < oo, then H is a tame isomorphism: —* ^{u}- 
If limsup = 0, then H is a tame isomorphism: <S( W ) — + A( W ). 



This theorem is proved in two parts, first the upper bound: 

Proposition 3.1. There exist C < oo, depending only on W, and C < oo (which may 
depend on uj) such that, for all f £ Sg u , 



||H(/)L<C||/Ilc w (9) 

and the lower bound: 



Proposition 3.2. Assume (v) There exist C < oo, depending only on w, and C < oo 
(which may depend on uj) such that, for all a € A^,, 



\\H- l (a)\l<C\\a\\c u . (10) 

The two inequalities above give immediately the projective case (jHJ, with |/| ■ := ||/L W and 
:= \\H(f)\\juj- The condition limsup = ensures that Proposition 13.21 can be 
applied to any jui. 

Replacing ui by C uj in the previous two propositions, we get 0, with |/| := \\f\\* 
and \\H(f)\\- := ||i?(/)|L. In that case, the condition limsup < oo guarantees that 
Proposition 13.21 can be applied to j when j is large enough. Thus Theorem [T] is proved. 

We now turn to the proof of Propositions 13.11 and 13.21 In the sequel, C denotes a constant 
(not necessarily persistent) which may depend on uj and C denotes a constant which may 
depend only on uj. In Lemmas 13.41 aud 13.51 K denotes a universal constant. 



10 

3.2. Upper bound 



JEAN-MARIE AUBRY 



The proof is based on a few elementary lemmata. We consider that a weight function uj has 
been fixed satisfying Definition [1] and we recall that uj := {Xuj, A > 0}. 

Lemma 3.3. Let f2 : t \— > w(e') and Jet fi* : v \— > sup t6R ^t — fl(t) be its convex conjugate. 
There exists a C < oo such that, for all f £ <S W , for al/ f > 0, 



and (here D = ■£:) 



\\x»fL><ctP*w\\f\\ 



\\D»f\\ L 2 < Ce n *^ +1) ||/|| 



Note that (ii) ensures f2*(i/) < oo for all f > 0. 



Proof. Since ||/|| w = ||/|| , wc only have to prove the first inequality. 
Note that f2* is increasing, so 



\x\ 2v \f{x)\ 2 dx<2sn V |/(.*)| 2 <2||/|| 2 

|as|<l |sc|<l 

< 2e- 2 ^( 1 'e 2n *^ +1) ||/||^. 



On the other hand wc have 



\x\ 2v \f(x)\ 2 dx < supe^+V^-^WfWl 

|x|>l x>l 

< 2 e 2sup *>o(' y+1 ) t ^ n O ||/|| 2 
<2e™*W\\f\C 



dx 



\x\>l 



Adding up the two, we get the result (with C = V2(e n ^> + 1)). 



□ 



The second lemma does the interpolation between the information on / and the information 
on /. 

Lemma 3.4. There exists a K < oo such that, for all weight functions uj, there exists C, 
for all f £ Su, for all li, v £ N, 

\\x v Wf\\v < CK»+"e = ll/IL- (11) 
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\x v D"f\\L' 



D tl f{x)x 2u D^f{x)dx 



f{x)D^(x 2 "D tl f(x))dx 



<£ 

k=Q 
k=0 



\xf\1v 



A (2i/)! 



fe=0 



A; / (2z/ — fc)! 

/LJI-d 2 "-"/ 



x 2u - k D 2 "- k f(x)dx 



k)\k 



j ^ e n* (2u~k+i)+n* (2/i-fc+i) || y || : 



thanks to Lemma 13731 and the fact that 6> w C 5. 

By convexity, v h-> - fi*(0) is superadditive, so fi*(2^- fc+l) < fi*(2^+ 1) - 



f2*(0), and the same with /x. Thanks to (iv) there exists xq =: e*° such that x > x$ =>• u>(x) < 
a; 2 . It follows that 

> sup At - Q(t) V sup fci - e 2t 

t«o t>to 

when k > ko '■= 2x\. In that case, using Stirling's formula, we see that there exists a constant 
B such that - n*(k) < fcf . 

At this point we have shown that when k > ko, 

log(fc!) + n*(2i/ -k+1) 

+Q*(2fi - k + 1) < Q*{2v + 1) + fi*(2^ + 1) + kB + 2fi*(0) 



using the superadditivity again, 

< VL*{2v + 2{jL + 2) + kB + 3Q*(0) 
wc crudely bound Of) ( 2 ^) by 4 M+I/ and e feB by (e 25 )^ 4 ^^, and the lemma follows. 



□ 



To finish, the proof of [2TJ Theorem 3.4] can be adapted, with Nl£ ■= (AQ)*(2(a + l))/2, 
to establish Proposition 13. II For the sake of completeness, we produce our own version of this 
proof. The third lemma uses the previous one to bound the iterates of the Hermite operator 
acting on /. 
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Lemma 3.5. There exists a K < oo such that, for all weight functions w, there exists C, 
for all f e S u , 

Proof. Using Leibnitz' rule, we can expand (x 2 — D 2 ) M as a sum of 2 M terms C^ M x v D^ 
with fj, + v < 2M, and each C M>V < 4 M . Applying (TTT]) then yields 

\\( X ^D 2 r f \\ L2 <C(SK) M e SJll ^\\f L 

which is the desired result. □ 



Proof of Proposition ^. 1\ Let u be fixed and let f £ S u . Recall that h n is also an eigenvector 
of the self-adjoint Hermitc operator x 2 — D 2 , with eigenvalue N := 2n + l. For M > 0, we thus 
have 



hence 



{(x 2 - D 2 ) M f, h n ) = </, (x 2 - D 2 ) M h n ) = N M (/, h n ) 



\{f,h n )\<N- M \\(x 2 -D 2 ) M f\\ L2 



using Lemma 



optimizing in A/ 



i sup M log(#)(2M+l)-n*(4M+2) iijii 



since f2 is convex 



< C (^Ve-K*^(*))||/| L 



finally using (ii) and (iii 



<Ce-^||/|| L 



for some C > that depends only on the constant in (iii) thus only on to. In the previous 
reasoning we can then replace ui by Coj to obtain □ 
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3.3. Lower bound 

To prove ([TO)) , we need a bound on the decrease of the Hermite functions. It is well known 
(see Szego's reference |25l 8.22.14] for instance) that |/i„(x)| < 1 for all x £ K, n £ N (actually 
\K(x)\ < n-T2). 

Lemma 3.6. Let y and z be two C 2 functions: [xq, +oo) — > (0, +oo) such that y'(x) — > 0. 
z is bounded, satisfying the differential equations 

y"{x) = (p(x)y(x) 
z"(x) = ip(x)z(x), 

with continuous <j)(x) < ip(x), and initial conditions y(xo) = z(xq). Then for all x > xq, 
z(x) < y(x). 

Proof. Suppose that there exists x\ > xq where z(x±) > y(x\). Then for some X2 £ [xq, x\] 
we have z'(x2) > y'{x2) and z(x2) > 2/(^2) ■ Consequently, for all x > X2, z"{x) — y"{x) > 0, 
and z'(x) — y'(x) > z'(x2) — y' (^2)- When x — > 00, liminf z'(x) > z'(x2) — y'ixz) > 0, which 
contradicts the boundedness of z. □ 

Lemma 3.7. For all n £ N and \x\ > s := \/2n + I, 

_ (|a;|-3) 2 _ (|x|-s) 2 

\h n (x)\ < h n (s)e 2 < e 2 . 

Proof. By parity, we can assume x > s (then h n (x) > 0). This implies that 

(x - s) 2 - 1< x 2 - s 2 . (12) 
Recall that h n satisfies the differential equation h!' n = (x 2 — s 2 )h n . On the other hand, 

(x-s) 2 

y{x) := h n (s)e 2 satisfies y" = ((x — s) 2 — l)y. This, the obvious properties of h n and y, 
and (Ti"2"|) together imply, by Lemma l3~6l that h n {x) < y{x). □ 

Proof of Proposition 13.21 Since limsup < |, there exists to < oo and 6* < 1 such that, 
for all t > t , u(t) < e 2t ^. Let n e (x) := 

We start with a £ (C to be determined later) and suppose that \\a\\ u < 1, which means 
that for all n > 0, |a„| < e-CMv 7 ™). Thus if / = ^3 a n h n (this series converging in iS), 

|/(x)|< e-^^\h n {x)\+ J2 e-^^\h n (x)\. 

0<n<ne(x) n>7io{x) 
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In the first sum, \x\ — y/2n + 1 > B\x\, so we can use Lemma 13.71 to bound |/i n (a;)| by 
e _e2 T" < e~ w ^ x ^ as soon as |x| > t . Because of the series Yln^ "^^™ 1 converges, so S± 
is bounded by Ce~ u ^ x \\ 

In the second sum, we simply bound |/i ra (x)| by 1 and use the following on the tail of the 
sum. 



n>y 



< 



-c^) dx 



e t-cn(t/2) dt 



< 



l°g(j/-l) 
1 



< 



C Q/( lo^) ) _ 1 
e log(y)-Cw( v ^) 



log(j/) 



l^t-cod)^ 



< 



(13) 



We apply this with y = n$(x), u>(y/y) > uj(^-^-\x\) > =co(\x\) by (iii)| for some C that depends 
only on ZJ. Finally we obtain Sa < Ce - "' ■ \ and the proposition is proved. 

□ 



Remark. Proposition 13.21 can also be shown by adapting the proof of I21[ Theorem 3.4] 



4. Wigner distribution 

We recall the definitions of & p and $(/, /), already given in the introduction as ([2]) and j3|. If 
p is a semi-definite positive Hermitian operator: L 2 — > L 2 that diagonalizes in an orthonormal 
basis ipi> with eigenvalues p^, then 



If pf represents a pure state, in other words, if it is the projector on the subspace generated 
by /, then we write 



2vr 



e^f( q -^)f(q+^)dx. 
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The squared modulus of a function / and its Fourier transform / can be recovered from 
<&(/, /) as marginal distributions 

|2 



and 



Let us define 



Hf,f)(q,p)dp=\f(q)\' 



m,f)(q,p)dq=\f(p)\ 



:=sup|$(q,p)|e 2(w(l « l)+ " (|p|)) 



(14) 



(15) 



and 



l*(/,/)L 

Clearly these are two norms, and using (|14[) and (I15|) we get by integration, for all q and for 
all p, 

\f(q)\e^ < V$*e- 2 " (M) MlfL 

as well as 

\f(p)\e^ < y/S^-^MfL 



the integrals being finite by (ii) hence finally 



< 



I 10 IIICJ 



(16) 



the constant C depending only on u>. We now aim at the corresponding lower bound. 



Proposition 4.1. Assume (v) There exists C < oo, depending only on ui, and C < oo 
(which may depend on u), such that for all f € 



Ill/Ill < c'WfWc,- 

We prove actually a more general result, which gives ([T7]) as a particular case. Let 

||HL:=sup|p m ,„|e^H-(v^). 

m.n 



(17) 



Proposition 4.2. Assume (v) There exists C < oo, depending only on u, and C < oo 
(which may depend on uj ), such that for all p, 



(18) 
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The proof is based on radial bounds for the special Hermite functions. We recall that 

m , n 

where $ miTl , defined by (JSJ), can also, as shown in |22| . be expressed as follows: when m > n, 

*m,n(?.P) = ^ (£]) (V2(ip - ofiMV + 2p 2 ). (19) 

Here := (n!) _1 e x x _Q ^r(e _;E a; n+Q ) is the (non normalized) Laguerre polynomial of degree 
n and order a. 

If to < n, then by Hcrmitian symmetry <fr m .n{q,p) = $„. m (g, — p), which is equivalent 
to taking the canonical generalization of Laguerre polynomials for — n < a < 0: '■= 
^-^- L (—x)~ a L~+ a (x). But since p is Hermitian, p n>m = p m , n , we only have to consider to > n 
in the sums below. 



The modulus of is thus radial: writing r := \J q 2 + p 2 , we have 

m—n 1 

l m ,n(r) ■= l*m,»(ff,p)l = — ( ~, V e~ r \ m ~ n \l% ~ n (2r 2 ) I . (20) 

7r \to! / 1 1 

What we need at this point is a bound on l m .n that is uniform on to and n, in the same 

fashion as Lemma [5771 



Lemma 4.3. There exists a constant K such that, for all m> n and s := y/m + n+l, for 
all r>0, 

!1 if0<r<s 
(21) 
e -( r - s ) 2 ifr>s. 

Proof. When r < s, the result follows from the uniform bounds on Laguerre polynomials, 
for instance given by Krasikov [19] : 

(L%(x)) 2 e- X x a+1 < 1444ri"5 (n + a + 1)* (22) 

to be used in ([20]) with x = 2r 2 and a = to — n. 

When r > s, L"(2r 2 ) doesn't vanish and keeps the same sign as L"(2s 2 ). Now, as it can 
be seen from [251 5.1.2], the function z(r) := y/rl min (r) satisfies the differential equation 
z " = (4(r 2 - s 2 ) + ° 2 ;, 1/4 )z. On the other hand, y(r) := ^l m , n {s)z~ (r ~ s)2 satisfies y" = 
(4(r - s) 2 - 2)y. When r > s. 

4(r- S ) 2 -2 <4(r 2 - S 2 ) + " 2 ~ 2 1/4 (23) 

from which we conclude with Lemma 13.61 that z(r) < y(r). □ 
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Proof of Proposition \4.1{ By (jU and Proposition l3.il if p = pt is a pure state, then for all 
m, n, 

\Prn,n\e^ )+ ^ = \ (f, h m ) |e<^ \ </, h n ) \ e"^ < C\\ f\$ u 

hence, applying Proposition ^. 2[ 



\HfJ)L<c 2 \\f\& 



C u' 

□ 



Proof of Proposition \4.2[ It is very similar to Proposition 13.21 There exists to < oo and 
9 < 1 such that, for all t > t , u(t) < 2 ^. Let m e {r) := r 2 (l — 0) 2 — 1 and assume that 
\\p\\cu < 1, C 7 to be chosen later. For all m,n we have |p roj „| < e -CMVS)+w(VH))_ When 
m + n< m e (r), r- s>6r and by ([21]), this means that l m>n {r) < Ke- 9 *^ < Ker 2 ^) . So 

J2 \Pm,n\lm,n( r ) < Ce~ 2 ^) ( 24 ) 
m+n<rag (r) 

for C := K J2 m n e^^^^f^+uiVn)) . 

On the other hand, comparing the sum to the integral, we get similarly as in (|13|) 

e -5Mv / S)+«(v / S)) < Ce""S w(v/I) 

m+n>y 



using (|2T|) again, 



m+n>me (r) 



< Ce- 2u){r) (25) 



if the constant C is chosed large enough (depending only on ui) . Combining (f2"4")l and (|25[) yields 
the announced result. 

□ 

If the exact form of ui is known, such as lo(x) — x 13 , < (3 < 2, a constant C close to the 
optimal can easily be obtained. An application to quantum statistics will be presented in a 
forthcoming paper. 

4.2. Mixed states 

If we apply (fl"6)) to the right-hand side of ([9]), then use ([4]), we obtain that the density matrix 
coefficients of a pure state / decrease as 

|p m , n |<C|||/|||^e-^^) + -(^) 
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in other words, we have a converse to Proposition ^. 21 but only when p = pf. No such converse 
can hold in the case of a general (mixed) state, as the following example shows. 
Following Butucea et al. [6j, let 



(-l) m S mn j(-ir S 1 z fn (l-z)dzi£m = n 
Pm ' n: -(m + D(m + 2) -| 0else _ 

We also recall an integral representation for Laguerre functions |251 5.4.1], for a > — 1: 



(26) 



ex ■ 



Then, using (fT9|) once again with a = m — n = and x := 2(g 2 +p 2 ), 
$ p (q,p) = — £ f z m (l - z)dzi° n (x) 



< 



7T 

e = 

7T 
7T 



V z m (l - z) — e- t J (2Vtx)dtdz 



(l-z) 



eS z - 1)t J {2Vtx)dtdz 



'.'-i-dz 



We see that, although p mi „ has algebraic decrease along the diagonal, & p (q,p) decreases like 

/ 2 , 2\ 

e -W +p )_ Cancellations occur because of the alternating signs in p m . m - 



4.3. Ambiguity function 

The ambiguity function A(f,f) of a signal f £ L 2 \s simply the Fourier transform in (q,p) 
of its Wigner transform. Knowing ([5]), 

- — (—j\ m + n 

A(f, /)(<?, 07) := *(/, /)(*, ttf) = £ ^ (/, M (/, M*m,n(f , f ) 

m,n 

The proof of Proposition 14. II works mutatis mutandis to obtain 



^4(/,/)L<qi/H- 



the constant C depending on uj and C depending on ui. 
Naturally, the equivalent of (Ti"T)|) holds also: 



<Cy/\A(f,f)\c u 



in particular <&(/, /) and /) are tamely equivalent in their respective ui norms. However, 
as the previous example (|26[) shows, this is not necessarily true for mixed states. 
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